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CHAPTER  I 
INTRODUCTION 

In  problems  of  water  impact,  the  compressive  stress  wave 
generated  in  a  weapon  vehicle  at  the  surface  of  contact  can  be 
quite  severe.  In  certain  instances  it  may  be  necessary  to  miti¬ 
gate  the  peak  stress  so  that  some  maximum  allowable  stress  is 
not  exceeded.  This  can  be  accomplished  by  protecting  the  contact 
surface  with  a  crushable  material.  Consider,  for  example,  the 
vertical  water  impact  of  a  torpedo  model  on  a  calm  water  surface. 

When  the  nose  enters  the  water  unprotected,  a  short-duration  high- 
amplitude  stress  pulse  is  generated  which  is  nearly  a  linear 
function  of  the  impact  velocity.  Figure  1.  However,  when  protected 
by  a  crushable  mitigator,  the  duration  of  the  stress  pulse  is 
increased  and  the  amplitude  is  decreased.  The  extent  of  mitigation 
can  be  substantial,  as  shown  in  Figure  1,  but  is  dependent  upon  the 
mechanical  properties  of  the  crushable  material  chosen.  For  a 
particular  weapon  design  and  impact  condition,  the  material  prop¬ 
erties  of  the  mitigator  should  be  optimized  on  the  basis  of  allowable 
load,  weight,  volume,  etc. 

Many  materials,  e.g.,  rigid  polyurethane  foam,  exhibit  elastic, 
nearly  perfectly  plastic  stress-strain  relations  until  compaction 
occurs  at  large  strains.  Compaction  refers  to  the  closure  of  the 
pores  in  the  foam  to  the  extent  that  it  behaves  as  a  solid  material. 
Stresses  greater  than  the  yield  stress  cannot  be  transmitted  by 
such  materials  until  compaction  occurs.  A  typical  compressive 
stress-strain  curve  for  rigid  polyurethane  foam  is  presented  in 
Figure  2 .  It  is  characterized  by  a  low  elastic  modulus  resulting 
in  a  gross  yield  strain  on  the  order  of  5  percent,  near  perfectly 
plastic  behavior  to  large  strain  on  the  order  of  50  percent,  and 
subsequent  increasing  rigidity  due  to  compaction.  The  unloading 
curve  is  characterized  by  an  initially  rapid  decrease  in  stress 
with  strain  followed  by  a  large  recovery  in  strain  at  very  low 
stress. 

The  purpose  of  this  report  is  to  investigate  the  mitigation  of 
shock  and  absorption  of  energy  under  impact- loading  conditions 
through  the  use  of  crushable  materials.  Experimental  data  are  pre¬ 
sented  on  the  mitigated  impacts  of  a  short  steel  projectile  on  a 
long  aluminum  rod.  A  theoretical  analysis  is  formulated  and  the 
results  are  compared  with  the  experimental  results.  Since  the  geom¬ 
etry  of  systems  subjected  to  impact  loading  is  often  dictated  by 
other  considerations,  mitigator s  of  varying  cross-sectional  area, 
such  as  conical  frustums,  as  well  as  mitigators  of  constant  cross- 
sectional  area,  are  included. 
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The  dynamic  properties  of  the  rigid  polyurethane  foam  were 
obtained  experimentally  and  used  as  input  for  the  theoretical 
analysis.  Functional  formulations  for  the  density  and  strain 
rate  dependence  of  the  elastic  modulus,  yield  strength  and  plateau 
stress  are  established  using  these  data  obtained  from  experiments 
and  from  the  literature.  These  empirical  functions  are  used  in  a 
computer  program  to  facilitate  the  prediction  of  mitigator  perform¬ 
ance  under  impact  loading. 

This  report  is  divided  into  six  sections: 

1.  In  Chapter  II,  the  theoretical  stress-strain  rela¬ 
tions  for  the  foam  are  formulated. 

2.  In  Chapter  III,  the  dynamic  material  property  data 
are  presented  and  the  theoretical  models  are  numerically  fitted 
to  these  data. 


3.  In  Chapter  IV,  the  analytical  treatment  of  the  impact 
problem  is  discussed. 

4.  In  Chapter  V,  the  apparatus  used  for  the  dynamic 
stress-strain  and  the  impact  experiments  are  described. 

5 .  In  Chapter  VI ,  the  theoretical  results  are  compared 
with  the  experimental  data. 

6.  In  Chapter  VII,  the  conclusions  are  presented. 


CHAPTER  II 
MATERIAL  MODELS 

Three  theoretical  stress-strain  models  are  considered,  Figure  3. 
The  first  is  simply  a  rigid,  perfectly  plastic,  locking  model.  The 
second  is  an  elastic,  perfectly  plastic,  compacting  model.  The  third, 
an  elastic,  perfectly  plastic,  strain-rate  dependent,  compacting 
model,  is  a  rate  dependent  version  of  the  second. 

RIGID,  PERFECTLY  PLASTIC,  LOCKING  MODEL 

The  rigid,  perfectly  plastic,  locking  model,  hereafter  referred 
to  as  the  RPL  model,  is  illustrated  in  Figure  3A.  It  ignores  all 
elastic  response  of  the  material  during  loading  and  unloading.  Any 
compressive  strain  up  to  the  compaction  strain  produces  a  stress 
equal  to  the  yield  stress  until  unloading  occurs.  Unloading  takes 
place  with  essentially  no  recovery  in  strain.  Perrone,  Reference  1, 
has  used  this  model,  with  the  addition  of  strain  rate  effects,  to 
obtain  a  closed  form  solution  for  the  dynamic  response  of  a  mass 
coupled  to  a  rigid  wall  by  a  wire.  Green,  et  al.,  Reference  2,  found 
the  yield  stress  for  rigid  polyurethane  foam  to  be  strain  rate  and 
density  dependent.  In  Reference  3,  Perrone  points  out  that  most  of 
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the  plastic  work  is  done  before  the  stress-strain  rate  point 
departs  appreciably  from  its  initial  value.  Therefore,  the 
yield  stress  associated  with  the  initial  strain  rate  can  be  used 
for  the  entire  calculation  with  little  loss  in  accuracy  and  much 
simplification. 

While  this  model  may  be  acceptable  for  predicting  plastic 
strain  and  energy  absorption  using  the  proper  initial  yield  stress, 
it  does  not  adequately  describe  the  force-time  and  velocity-time 
curves  generated  during  impact. 


ELASTIC,  PERFECTLY  PLASTIC,  COMPACTING  MODEL 


The  elastic,  perfectly  plastic,  compacting  model,  the  EPC 
model,  includes  the  addition  of  pre-yield  elastic  response,  elastic 
unloading  and  post-compaction  elastic  response,  Figure  3B.  Schmidt 
and  Hoffman,  Reference  4,  and  Mustin,  Reference  5,  have  presented 
the  rigidity  modulus  for  polyurethane  foams  as  a  function  of  tem¬ 
perature,  Figure  4.  These  data  indicate  that  since  the  rigid  foams 
are  in  the  glassy  state  at  room  temperature,  the  rigidity  modulus 
is  nearly  temperature  independent  at  that  temperature.  Schwaber 
and  Meinecke,  Reference  6,  suggested  that  since  these  foams  are  in 
the  glassy  state,  the  modulus  is  also  approximately  strain  rate 
independent.  However,  the  modulus  is  highly  density  dependent  and 
has  been  extensively  studied  under  "quasi-static"  loading  by  many 
of  the  manufacturers.  Schwaber  and  Meinecke,  Reference  6,  also 
stated  that  the  unloading  stress-strain  curve  appears  to  be  inde¬ 
pendent  of  strain  rate.  The  slope  of  the  unloading  path  is  assumed 
to  be  equal  to  the  slope  of  the  initial  elastic  loading  path  because 
reliable  experimental  data  on  the  actual  unloading  path  are  lacking. 
The  large  recovery  in  strain  at  very  low  stress.  Figure  1,  is  ignored 
also  because  of  insufficient  data.  However,  this  is  not  a  serious 
deficiency  except  for  multiple  impacts.  No  multiple  impacts  are 
considered  in  this  work.  The  post-compaction  response  is  treated 
elastically  here  for  convenience.  In  reality,  the  compaction  process 
does  not  occur  abruptly  and  the  compaction  strain  is  therefore  a 
nebulous  quantity,  Figure  2.  The  transition  is  smooth  and  amenable 
to  modeling,  but  data  are  not  available  on  the  density  and  strain 
rate  dependence  of  the  process.  In  general,  the  loading  and  unloading 
paths  following  the  onset  of  compaction  are  not  coincident. 


ELASTIC t  PERFECTLY  PLASTIC,  STRAIN  RATE  DEPENDENT,  COMPACTING  MODEL 

The  elastic,  perfectly  plastic,  strain  rate  dependent,  com¬ 
pacting  model,  the  EPRC  model,  introduces  a  strain  rate  dependent 
plateau  stress  in  addition  to  a  strain  rate  dependent  yield  stress, 
Figure  3C.  The  elastic  response  remains  strain  rate  independent, 
while  the  yield  stress  and  yield  strain  become  strain  rate  dependent, 
increasing  with  increasing  rate.  The  plateau  stress  becomes  a 
function  of  the  instantaneous  strain  and  strain  rate.  It  will 
therefore  differ  from  the  yield  stress  for  situations  of  varying 
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strain  rate.  The  unloading  path  is  again  chosen  as  elastic  with 
a  slope  equal  to  that  of  the  initial  elastic  loading.  However, 
for  cases  of  equal  maximum  strain,  the  path  is  not  unique  except 
for  cases  of  equal  strain  rate.  Again,  for  convenience,  the  com¬ 
paction  strain  is  chosen  as  a  discrete  strain  regardless  of  strain 
rate.  The  post-compaction  loading  and  unloading  path  is  also  not 
unique  except  for  cases  of  equal  strain  rate. 


CHAPTER  III 

NUMERICAL  FIT  OF  MODELS  TO  DATA 

For  the  analysis,  the  following  assumptions  are  made  concerning 
the  polyurethane  foam  material. 

1.  The  elastic  modulus  is  density  dependent,  but  strain  rate 
independent . 

2.  The  yield  stress  and  plateau  stress  are  density  dependent 
and  strain  rate  dependent. 

3.  The  unloading  stress-strain  relation  is  density  dependent, 
but  strain  rate  independent. 

Material  property  data  for  rigid  polyurethane  foam  have  been 
obtained  from  stress-strain  strain-rate  tests  and  from  the  literature. 
Since  foam  materials  are  typically  anisotropic  at  low  densities, 
all  data  and  test  results  apply  to  material  loaded  parallel  to  the 
direction  of  rise.  The  direction  of  rise  is  defined  as  the  direction 
in  which  the  foam  is  allowed  to  expand  during  fabrication. 

As  previously  stated,  the  elastic  modulus  is  assumed  to  be 
density  dependent,  but  strain  rate  independent.  Therefore,  data 
acquired  under  "quasi-static"  conditions  are  considered  applicable 
to  strain  rates  up  to  at  least  1000  in. /in. /sec.  Modulus  data 
obtained  according  to  ASTM  D1621  are  presented  in  Reference  7  as  a 
function  of  density.  These  data  are  listed  in  Table  1.  To  facilitate 
the  numerical  computations, a  curve  of  the  form 

_  _  k 

E  =  CP  (1) 

was  fitted  to  these  data.  The  resulting  curve  fit,  Figure  5,  is 

E  =  230  p1  * 7  (2) 

where  E  is  the  elastic  modulus  in  psi  and  p  is  the  initial  foam 
density  in  lbs./ft.^. 

The  dynamic  yield  stress  and  the  plateau  stress  at  20  percent 
strain,  Figure  2,  were  determined  experimentally  as  functions  of 
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TABLE  1 

COMPRESSIVE  MODULUS  OF  RIGID  POLYURETHANE  FOAM 

Tested  According  to  ASTM  D1621  at  74°F 
in  Direction  Parallel  to  Rise 


Density  Compressive  Modulus 

(lbs. /ft.3)  (psi) 


3 

4 
6 
8 

10 

12 

14 

16 

20 


1,020 

2,100 

4,400 

7,700 

11,500 

16,000 

21,000 

26,000 

38,000 


From  Reference  7 
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density  and  strain  rate.  Three  densities  of  the  foam  -  nominal 
3.5,  7.0,  and  14.5  lbs. /ft.  -  were  tested  in  compression  at  rates 
of  strain  between  10”^  and  10^  in. /in. /sec.  All  of  the  test  speci¬ 
mens  were  loaded  parallel  to  the  direction  of  rise.  Each  test  was 
conducted  at  nearly  constant  strain  rate.  Varying  densities, 
nonuniform  cellular  structures  and  flaws  produced  scatter  in  the 
data.  Typical  stress-strain  curves  are  presented  for  various  strain 
rates  in  Figures  6,  7,  and  8  for  each  of  the  three  densities.  Note 
the  drop  in  stress  which  occurs  after  yield.  Since  the  material 
models  omit  the  peak,  a  plateau  stress  -  for  example,  the  stress  at 
20  percent  strain  -  is  a  better  parameter  to  use  in  the  analysis 
than  the  peak  stress.  To  aid  in  determining  the  resultant  trends, 
equations  of  the  form 


o  = 


(cx  +  c2  In 


(3) 


were  fitted  to  the  yield  stress  and  plateau  stress  data  by  the 
method  of  least  squares.  £  is  the  average  strain  rate  in  the 
specimen  in  in. /in. /sec. ,  p  is  the  initial  density  of  the  specimen 
in  lbs. /ft. ^  and  C^,  C 2,  and  kn  are  the  parameters  adjusted  to 
satisfy  the  least  square  condition.  The  equation  developed  to  fit 
the  yield  stress  data  is 


oy  =  (10.9  +  .385  In  e)p1,64  (4) 

and  the  equation  developed  to  fit  the  plateau  stress  data  is 

op  =  (8.69  +  .306  In  e)p1,69  (5) 


The  density  and  strain  rate  dependence  of  these  data  and  the  fitted 
curves  are  illustrated  in  Figures  9,  10,  11,  and  12. 


CHAPTER  IV 
ANALYTICAL  APPROACH 

In  formulating  the  analytical  approach  to  this  impact  problem, 
the  following  assumptions  are  made: 

1.  The  projectile  is  sufficiently  rigid/  compared  to  the 
mitigator, that  it  can  be  treated  as  a  lumped  mass  which  undergoes 
negligible  deformation. 

2.  The  mass  of  the  mitigator  and,  hence,  the  inertia  is 
negligible  compared  to  that  of  the  projectile. 
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Assumption  (2)  implies  that  the  theoretical  results  will  be  the 
same  for  the  case  of  a  mitigator  attached  to  the  projectile  and 
the  case  of  a  mitigator  attached  to  the  rod. 


CONSTANT  AREA  MITIGATOR 


The  experimental  configuration  can  be  divided  into  three  free 
body  diagrams. 


- >+  u 

where 


o 


Summing  the  forces  applied  to  the  projectile  in  the  axial  or 
x-direction  gives 


-F  =  Ma 


(6) 


or 


(7) 


where ,  in  general , 
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F  =  F(e,  e,  A) 

e  *  strain  of  the  mitigator  (8) 

e  =  strain  rate  of  the  mitigator 

A  =  stress  area  of  the  mitigator 

Integrating  with  respect  to  time  gives  the  velocity  relation 

t 


o 


Integrating  again  gives  the  distance  traveled 


X1  *  uoT  - 


II  * 


it  dtdt 


by  the  projectile. 


Utilizing  the  one-dimensiona^  theory  of  finite  amplitude  wave 
propagation,  the  transient  response  of  the  rod  to  a  force  F  applied 
to  the  end  is  given  by 


Integrating  yields 


(12) 


.The  solution  to  this  impact  problem  lies  in  determining 
F(e,e,A)  at  each  instant  in  time.  To  obtain  this  result,  the  strain 
and  strain  rate  of  the  mitigator  must  also  be  determined  at  each 
instant  in  time.  The  engineering  strain  of  the  mitigator  is  given  by 
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=  L  -  (X2  -  Xx) 


(13) 


Substituting  for  X^  and  X2  yields 

T  t 


e  =  I 


u° '-//  ■  -«*-/  F 


o  o 


tugsr— pt 


(14) 


The  strain  rate  of  the  mitigator  is  given  by 


e  =  -<U2-U1> 

L 


(15) 


Substituting  for  and  U2  yields 

r  t 


o 


(16) 


Using  Equation  (8)  as  a  constraint,  Equations  (14)  and  (16)  can  then 
be  integrated,  see  Appendix  A,  to  determine  the  force-time  curve,  the 
strain-time  curve  of  the  mitigator,  and  the  velocity-time  curves  of 
the  projectile  and  the  rod. 


Of  equal  interest  to  the  magnitude  of  the  force  transmitted  by 
the  mitigator  is  the  amount  of  energy  absorbed  by  the  mitigator. 

For  example,  a  soft  elastic  material  will  not  dissipate  as  much 
energy  as  a  soft  plastic  material.  The  strain  energy  absorbed  by 
the  mitigator,  plus  the  kinetic  energy  and  the  strain  energy  imparted 
to  the  rod,  equals  the  change  in  kinetic  energy  of  the  projectile 
during  the  impact.  The  final  strain  energy  absorbed  by  the  mitigator 
is  then 


where 


E.  =  ARE  .  -  ARE  -  ASE 
A  pro}  rod  rod 


2  2 

ARE  .  =  4  M(U  -  Ur  ) 
proj  2o  r 


(17) 


(18) 


and  from  Reference  8 


rod  =  2  <M*>rod  f  u2<it 

*o 


(19) 
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For  a  semi-infinite  rod  stressed  below  the  elastic  limit 


E2e2 


The  change  in  kinetic  energy  becomes 


I  (AEgl<*  f 


e  j  dt 
rod 


The  strain  energy  imparted  to  the  rod  is 


■if 


aedV 


Substituting  E  e  for  a  and  Agdt  for  dV  gives 


ASErod  ‘  5  <“9>roa  |  croddt 


Equation  17  then  becomes 


EA  -  I  «(°I  -  ) 


/  £r=a  dt 


ea  -  I  M(uJ  -  oj) 


<a )  f 

AE'rod  j 

~r\ 


F2  dt 


VARYING  AREA  MITIGATOR 

For  many  applications  the  mitigators  are  not  of  constant  area, 
but  rather  are  elements  of  varying  area.  Consider  the  conical 
frustum  in  Figure  13  to  which  an  increasing  load  F(t)  is  applied. 

It  is  divided  into  N  cylindrical  discs  of  increasing  cross-sectional 
area.  Now  consider  just  the  disc  N  =  1.  It  undergoes  elastic 
deformation  until  it  reaches  its  yield  strength,  then  perfectly 
plastic  deformation  until  compaction  occurs.  At  the  same  time, 
discs  N  =  2  through  N  =  N  undergo  only  elastic  deformation  because 
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of  their  larger  stress  areas.  Following  compaction,  disc  N  =  1 
can  transmit  a  higher  load.  Disc  N  =  2  then  yields,  deforms 
plastically  and  compacts.  This  process  continues  until  either  the 
load  is  restricted  or  all  of  the  elements  have,  in  turn,  yielded 
and  compacted. 

For  stresses  below  the  yield  stress  of  disc  N  =  1,  the  load¬ 
carrying  capability  of  the  mitigator  is 

F  =  OjAi  (26) 

where  and  A..  is  the  cross-sectional  area  of  disc  N  =  1.  For 

stresses  above  the  yield  stress  of  disc  N  =  1 ,  but  below  that  of 
disc  N  =  N 


F  =  aYh±  (27) 

where  is  the  cross-sectional  area  of  the  disc,  which  is  just 
loaded  to  its  yield  stress  with  limits 


The  area  A^  is  given  by 


Ai  =  J  D? 


where  is  the  diameter  of  the  mitigator  at 


(28) 


Xi 


(29) 


and  and  ec  are  defined  as  the  average  mitigator  strain  and 

the  compaction  strain  of  the  material,  respectively.  For  conical 
mitigators 


Dx  +  (D2  -  DX) 


(30) 


or 


Di  =  D-l  +  (D2 


The  area  A.  then  becomes 

l 


A-:  = 


ir 

? 


C1  +  ^2  ”  Dl^ 


I 


avgl 


(31) 


(32) 
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.A  compaction  strain  ec 
the  rigid  polyurethane 
Equation  (32)  gives 


of  0.50  in. /in.  appears  reasonable  for 
foam.  Inserting  this  information  into 

r  *12 


T 


D1  +  (d2 


(33) 


For  stresses  above  the  yield  stress  of  disc  N  =  N 


F 


0..A 
N  N 


(34) 


Equations  (14)  and  (16)  can  now  be  integrated  in  exactly  the  same 
manner  given  for  the  constant  area  case  using  Equation  (26),  (27), 
or  (34)  to  determine  the  appropriate  force  F. 


CHAPTER  V 

EXPERIMENTAL  APPARATUS 
IMPACT  TESTING  APPARATUS 


The  impact  experiments  were  conducted  using  the  apparatus 
illustrated  schematically  in  Figures  14  and  15.  Short  -  1/2,  3/4  and 

I- inch-long  -  cylindrical  steel  projectiles  were  accelerated  to 
the  impact  velocity  using  a  low-pressure  air  gun.  The  target,  an 

II- foot-long  aluminum  rod  supported  on  Teflon  mounts  to  allow  free 
axial  movement,  was  placed  on  axis  with  the  gun  barrel  at  the  muzzle. 
The  projectiles  and  the  rod  were  nominally  1/2  inch  in  diameter. 

The  mitigators  were  bonded  to  the  impact  end  of  the  rod  with  a 
thin  film  of  grease. 

In  order  to  obtain  reasonable  planarity  between  the  impact 
faces  of  the  projectile  and  the  mitigator,  the*gun  bore  and  the 
rod  were  aligned  with  the  aid  of  a  telescope.  In  addition,  the 
impact  was  arranged  to  occur  while  the  projectile  was  partially 
within  the  gun  barrel.  The  driver  gas  was  vented  through  holes  in 
the  barrel  to  eliminate  the  driving  force  from  the  projectile 
before  the  impact  occurred.  Thus,  free-f light  conditions  existed 
in  the  axial  direction  at  the  time  of  impact. 

Two  light-sensitive  transducers  were  placed  one  foot  apart 
along  the  gun  barrel  near  the  muzzle  to  record  the  projectile 
transit  time  over  that  interval  before  impact.  A  small  hole  was 
drilled  in  the  barrel  at  each  location  to  allow  light  to  impinge 
upon  the  transducer.  As  the  projectile  passed  each  location  and 
the  light  beam  was  interrupted,  each  transducer  generated  an  elec¬ 
trical  signal.  These  signals  were  used  to  start  and  stop  a  Beckman/ 
Berkeley  Model  7260R  time  interval  counter.  The  impact  velocity 
was  then  determined  from  the  relation 
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U  _  AX  _  12  (inches)  ,35j 

°  AT  -  AT  (seconds) 

A  Hycam  high-speed  framing  camera  (5000  f rs\mes/sec . )  was  used 
to  record  the  maximum  strain  of  the  mitigators.  Although  the  entire 
time  of  contact  between  the  projectile  and  the  mitigator  typically 
covers  only  5  to  10  frames,  the  accuracy  of  the  technique  is  reason¬ 
ably  good.  Since  the  rate  of  change  in  strain  is  low  at  the  time 
of  maximum  strain,  the  largest  error  encountered  in  the  maximum 
strain  measurement  is  less  than  3  percent.  The  projectile  rebound 
velocity  was  also  determined  from  these  films  from  the  relation 

Ur  =  (distance/frame)  .  (frames/second)  (36) 

Strain  gages  were  installed  on  the  aluminum  rod  12  inches 
from  the  impact  end  to  record  the  force  pulse  transmitted  by  the 
mitigator.  The  gages  were  so  installed  and  connected  in  a  Wheatstone 
bridge  as  to  cancel  out  the  effect  of  bending.  An  excitation 
voltage  of  6  volts  was  used  to  power  the  bridge.  Prior  to  each  test 
the  bridge  was  balanced.  During  each  test,  the  bridge  imbalance 
was  recorded  on  a  Textronix-type  531  oscilloscope,  fitted  with  a 
type-D  plug-in  unit  and  a  type  123  preamplifier.  The  gain  of  the 
preamplifier  was  100  to  1.  The  oscilloscope  traces  were  recorded 
on  Polaroid  film  and  reduced  to  digital  form  using  a  telereader. 

For  these  particular  test  conditions,  the  11-foot  rod  was  of  suffi¬ 
cient  length  to  allow  the  transmitted  pulse  to  completely  pass  the 
gage  location  before  the  leading  edge  of  the  pulse  reflected  from 
the  free  end  of  the  rod  arrived  at  the  gage  location.  Hence,  no 
analysis  of  the  reflected  pulse  was  necessary  to  unfold  the  trans¬ 
mitted  pulse  from  the  data. 


COMPRESSION-TESTING  APPARATUS 


As  input  for  the  theoretical  analysis,  "quasi-static,"  medium 
strain  rate  and  high -strain -rate  stress-strain  curves  for  the  rigid 
polyurethane  foam  were  experimentally  determined  in  compression  using 
a  Baldwin-Southwark  universal  testing  machine,  a  Plastechon  Model  581 
medium-strain-rate  testing  machine  and  a  split  Hopkinson  pressure 
bar,  respectively.  The  Baldwin-Southwark  machine  was  a  standard 
hydraulic-type  machine  with  graphical  load-deflection  readout.  The 
Plastechon  Model  581  was  simply  a  pneumatic-hydraulic  version  of 
the  standard-type  machine  capable  of  higher  crosshead  speeds,  but 
without  servo  control.  This  machine  was  not  sufficiently  rigid  to 
measure  strain  via  measuring  crosshead  displacement.  Therefore,  a 
pair  of  matched  Optron  Model  680AX  electro-optical  trackers  were 
used  to  measure  the  strain  by  tracking  the  displacements  of  the  two 
faces  of  the  specimens.  The  tracker  outputs  were  recorded  on  an 
oscilloscope  equipped  with  a  differential  amplifier.  The  load  was 
simultaneously  recorded  on  the  second  channel  of  the  same  oscillo¬ 
scope.  The  oscilloscope  traces  were  photographically  recorded, 


13 


NOLTR  73-162 


converted  to  digital  form  with  the  aid  of  a  telereader  and  combined 
to  form  the  stress-strain  curves.  The  corresponding  strain  rate  was 
determined  from  the  slope  of  the  strain- time  data. 

The  split  Hopkinson  pressure  bar  Wc?  used  to  obtain  stress- 
strain  data  at  strain  rates  on  the  order  of  1000  in. /in. /sec.  The 
basic  concepts  which  govern  the  use  of  the  split  Hopkinson  pressure 
bar  are  found  in  References  9  and  10.  The  material  specimen  is 
placed  between  two  elastic  bars.  A  compressive  stress  wave  is 
propagated  along  one  bar  to  the  bar  specimen  interface.  This  stress 
wave  is  partially  reflected  and  partially  transmitted  through  the 
specimen  into  the  second  bar.  Taking  compressive  strains  as  positive, 
the  average  strain  rate  in  the  specimen  is  related  to  the  reflected 
strain  in  the  first  bar  as  follows: 


(37) 


where  g  is  the  wave  velocity  in  the  bar  and  L  is  the  initial  length 
of  the  specimen.  The  specimen  strain  may  be  determined  by  the 
integration  of  the  above  relation. 


«s  -"§*/  er  dt  W> 

O 

The  average  specimen  stress  is  related  to  the  transmitted  strain 
as  follows: 


o 


s 


E 


(39) 


where  E  is  the  Young's  modulus  of  the  bar  material  and  A^/As  is  the 
area  ratio  of  the  bar  to  the  specimen. 


The  split  Hopkinson  pressure  bar  experiments  were  performed 
at  the  same  facility  as  the  impact  experiments.  The  long  rod  was 
replaced  by  the  two  steel  pressure  bars.  The  short  projectiles 
were  replaced  by  a  10-inch-long  projectile  to  produce  an  appropriate 
stress  wave. 


CHAPTER  VI 

EXPERIMENTAL  AND  THEORETICAL  RESULTS 

Prior  to  conducting  each  experiment,  the  dimensions,  mass  and 
density  of  the  mitigator  to  be  tested  were  determined.  During  each 
experiment,  the  following  quantities  were  measured: 
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1.  Projectile  impact  velocity 

2.  Projectile  rebound  velocity 

3.  Maximum  mitigator  strain 

4.  Transmitted  force  as  function  of  time 

Following  each  experiment,  the  initial  projectile  energy  and  the 
energy  absorbed  by  the  mitigator  were  calculated.  These  data  are 
presented  in  Table  2. 

Typical  data  records  of  the  strain  in  the  rod  are  shown  in 
Figure  16  for  nominally  18.5  lbs. /ft.  cylindrical  and  conical 
foam  mitigators  impacted  at  approximately  410  in. /sec.  by  a  0.05- 
pound  projectile.  Since  the  rod  remains  elastic  at  all  times,  the 
strain  in  the  rod  is  directly  proportional  to  the  force  transmitted 
by  the  mitigator.  Therefore,  these  strain-time  data  can  be  viewed 
as  force-time  data.  The  curves  are  characterized  by  steep,  but 
finite  rises,  slowly  decreasing  or  increasing  plateaus, and  gradual 
decays.  For  the  cylindrical  mitigator,  the  force  decreases  slightly 
along  the  plateau  due  to  the  strain  rate  effect.  For  the  conical 
mitigator,  the  effect  due  to  the  area  change  is  sufficient  to  offset 
the  effect  due  to  strain  rate.  The  continual  increase  in  stress 
area  causes  the  force  to  increase  slightly  along  the  plateau. 

The  critical  test  of  the  model  is  its  ability  to  match  the 
transmitted  force  pulse  measured  in  the  experiments.  Energy  and 
velocity  comparisons  are  less  sensitive  to  the  details  of  the 
modeling.  Figure  17  illustrates  the  application  of  the  three  models 
to  the  case  of  the  0.05-pound  projectile  impacting  a  3/8-inch-diameter 
by  3/8-inch-long,  18.5  lbs. /ft. 3  foam  cylindrical  mitigator  at 
410  in. /sec.  Clearly,  the  RPL  model  does  not  predict  a  realistic 
force  pulse.  The  rise  time  is  infinitesimal,  the  plateau  is  flat, 
and  the  delay  time  is  infinitesimal.  The  addition  of  elasticity  in 
the  EPC  model  introduces  a  finite  rise  time  and  somewhat  longer 
decay  time  in  the  calculated  force  pulse.  Except  for  the  flat 
plateau,  this  force  pulse  resembles  the  experimental  data  very 
closely.  Adding  the  strain-rate-dependent  plateau  stress  in  the 
EPRC  model  produces  the  decreasing  force  along  the  plateau  seen  in 
the  experimental  data. 

Similarly,  Figure  18  illustrates  the  application  of  the  three 
models  to  the  case  of  the  0.05-pound  projectile  impacting  a  3/8-inch- 
diameter  by  1/2-inch-diameter  by  3/8-inch-long,  18.5  lbs. /ft.  foam, 
conical  mitigator  at  410  in. /sec.  The  results  are  similar  to  those 
obtained  for  the  cylindrical  mitigator  except  for  the  increasing 
force  along  the  plateau.  It  is  interesting  to  note  that  because  of 
the  competing  area  and  strain  rate  effects,  the  EPRC  model  does  not 
predict  the  peak  force  to  occur  at  maximum  strain.  It  will  be  seen 
later  that  this  result  agrees  well  with  the  experimental  data. 
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Although  the  EPRC  model  produces  the  most  representative 
transmitted  force-time  curve,  a  dip  develops  in  the  calculated 
curve  as  the  projectile  comes  to  rest  and  begins  to  rebound.  The 
dip  is  not  evident  in  the  experimental  data.  At  this  point  in  the 
calculation,  the  strain  rate  is  decreasing  very  rapidly  and  the 
predicted  strain-rate  dependence  of  the  plateau  stress  is  excessive. 
Since  the  dynamic  material  property  data  were  obtained  under 
constant  strain-rate  conditions,  and  constant  strain-rate  conditions 
do  not  exist  under  impact  loading,  these  data  are  not  expected  to 
apply  in  every  detail.  More  information  is  needed  on  the  time 
dependence  of  the  mechanical  properties  of  foam  materials. 

Furthermore,  the  logarithmic  function  fitted  to  the  stress- 
strain  rate  density  data  (Equation  (3)) 

.  K 

a  =  (C^  +  C2  In  e) p 


becomes  undefined  at  zero  strain  rate.  The  problem  is  averted  in 
the  calculation  simply  by  allowing  the  time  step  to  bridge  the  zero 
strain-rate  point.  Other  investigators,  for  example,  Perrone,  have 
used  a  different  function 


to  describe  the  strain-rate  dependence  without  encountering  the 
above  problem.  Instead  the  problem  becomes  that  of  determining 
a0.  As  applied  to  foam  materials,  the  authors  do  not  believe 
that  oQ  is  statically  determinable.  To  avoid  the  aforementioned 
dip  in°the  force-time  curve,  the  value  of  oQ  must  be  greater  than 
that  determined  statically.  Unless  the  proper  dynamic  0 p  is 
obtained,  Equation  (40)  offers  no  improvement  over  Equation  (3) . 

Theoretical  results,  corresponding  to  the  experiments,  were 
obtained  using  the  computer  program  listed  in  Appendix  B.  Impact 
velocity,  projectile  mass,  foam  density  and  mitigator  length  were 
varied  to  test  the  flexibility  of  the  EPRC  model.  Maximum  miti¬ 
gator  strain  and  energy  absorption  are  given  in  Table  2  for  comparison. 

The  experimental  and  theoretical  transmitted  force-time  curves 
of  test  1  are  presented  in  Figure  19.  These  results  are  for  an 
18.3  lbs. /ft. ^  cylindrical  mitigator.  The  experimental  curve  shows 
considerably  more  structure  than  the  theoretical  analog.  The  initial 
peak  may  be  caused  by  the  tendency  of  this  material  to  peak  at  yield 
and  then  drop  rapidly  to  some  plateau  stress  at  the  high  strain  rate. 
Subsequent  oscillations  are  due  more  to  electronic  noise  than  to 
material  response.  In  general,  the  theory  matches  the  experiment 
well  with  slight  differences  in  amplitude  and  pulse  length. 
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Similar  curves  for  test  2  are  shown  in  Figure  20.  These 
results  are  for  an  18.9  lbs. /ft. 3  conical  mitigator  subjected  to 
nearly  the  same  impact  conditions  as  those  of  test  1.  The  agree¬ 
ment  between  experiment  and  theory  is  excellent/  except  for  a  small 
difference  in  amplitude.  Note  that  the  peak  in  the  theoretical 
curve,  which  occurs  before  the  mitigator  reaches  maximum  strain, 
parallels  the  experimental  data  quite  well. 

Tests  similar  to  tests  1  and  2  were  performed  on  nominally 
10.8  lbs. /ft. 3  foam  mitigators  to  determine  the  effect  of  miti¬ 
gator  density  or.  the  transmitted  force-time  curve  and  energy 
absorption.  The  curves  for  tests  3  and  4  are  given  in  Figures  '21 
and  22,  respectively.  Somewhat  poorer  agreement  in  amplitude 
between  theory  and  experiment  is  encountered.  The  magnitude  of 
the  transmitted  force  is  about  40  percent  of  that  transmitted  by 
an  18.6  lbs. /ft.  foam  mitigator.  The  pulse  length,  on  the  other 
hand,  is  more  than  doubled. 

Tests  5  and  6  were  performed  under  the  same  impact  conditions. 

The  mitigators  used  in  these  tests  were  of  the  same  nominal  density 
as  those  of  tests  3  and  4,  but  were  30  percent  longer.  The  trans¬ 
mitted  force-time  curves  for  these  tests  are  shown  in  Figures  23 
and  24.  The  results  are  not  significantly  different, except  for 
smaller  maximum  strains,  Table  2. 

For  tests  7  and  8  the  impact  velocity  was  reduced  to  210  ft. /sec. 
All  other  aspects  of  these  tests  remained  similar  to  those  of  tests 
5  and  6.  The  resulting  force-time  curves  are  presented  in  Figures  25 
and  26.  The  rise  times  of  the  pulses  have  been  lengthened  and  the 
discrepancy  between  the  experimental  and  theoretical  results  has 
increased.  The  pulse  amplitudes  are  lower  because  of  strain-rate 
effects  and  compare  more  favorably  than  those  of  tests  5  and  6.  The 
pulse  lengths  are  shorter  and  agreement  is  good  between  theory  and 
experiment. 

Returning  to  the  use  of  a  short  mitigator,  test  9  was  performed 
in  a  similar  fashion  to  test  8  with  similar  results,  Figure  27.  The 
theory  predicts  an  increase  in  amplitude  due  to  strain-rate  effects 
with  a  corresponding  reduction  in  the  pulse  length.  However,  the 
experimental  data  from  tests  8  and  9  are  virtually  the  same.  Varia¬ 
tions  in  material  properties  due  to  nonhomogeneity  are  probably  suf¬ 
ficient  to  overshadow  the  predicted  differences  between  these  tests. 

Tests  10  and  11  are  similar  to  test  4  with  the  exception  of  a 
change  in  the  projectile  weight.  For  test  10,  the  weight  was 
reduced  by  25  percent  and  for  test  11,  the  weight  was  reduced  by 
50  percent.  Since  the  impact  velocity  is  essentially  the  same,  the 
kinetic  energies  of  the  projectiles  are  proportional  to  their  relative 
weights.  And,  since  the  force  developed  upon  impact  is  essentially 
the  same,  the  major  effect  of  reducing  the  projectile  weight  is  to 
shorten  the  pulse  length,  Figures  28  and  29. 
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The  experimental  force-time  curve  for  test  12  is  presented  in 
Figure  30.  The  impact  velocity  for  this  test  was  increased  to 
618  in. /sec.  For  this  projectile  velocity,  the  mitigator  should 
have  compacted.  However,  the  foam  material  crumbled  or  exploded 
due  to  the  air  trapped  within  the  voids.  This  explains  why  no 
large  increase  in  force  due  to  compaction  is  seen  in  the  data.  At 
low  strain  rates,  the  foam  would  have  simply  crushed.  Apparently, 
the  maximum  amount  of  energy  that  a  rigid  foam  mitigator  can  absorb 
without  self-destructing  is  dependent  upon  strain  rate,  decreasing 
with  increasing  rate. 


CHAPTER  VII 
CONCLUSIONS 

The  experimental  data  and  the  theoretical  results  contained 
in  this  report  permit  the  following  conclusions  regarding  the 
validity  of  the  material  model  and  the  effectiveness  of  rigid 
polyurethane  foam  as  a  shock  mitigator. 

1.  The  elastic,  perfectly  plastic,  strain-rate  dependent 
compacting  model  adequately  predicts  the  force  transmitted  by  the 
mitigator  as  a  function  of  time  and  the  energy  absorbed  by  the 
mitigator. 


2.  The  elastic,  perfectly  plastic,  strain-rate  dependent 
compacting  model  is  an  improvement  over  the  commonly  used  rigid, 
perfectly  plastic  locking  model. 

3.  Mitigators  of  varying  cross-sectional  area  are 
appropriately  modeled. 

4 .  More  information  is  needed  on  the  time  dependency  of 
the  plateau  stress  of  foams. 

5.  Rigid  polyurethane  foam  is  an  effective  shock-miti¬ 
gating  material,  Figure  1,  as  long  as  compaction  does  not  occur. 
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FIG.  2  TYPICAL  STRESS- STRAIN  CURVE  FOR  RIGID  POLYURETHANE  FOAM 
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A-RIGID,  PERFECTLY  PLASTIC,  LOCKING  MODEL 


B-ELASTIC,  PERFECTLY  PLASTIC,  COMPACTING  MODEL 


C-ELASTIC,  PERFECTLY  PLASTIC,  STRAIN  RATE  DEPENDENT,  COMPACTING  MODEL 


FIG.  3  MATERIAL  MODELS 
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TEMPERATURE  (°C) 

FIG.  4  RIGIDITY  MODULUS  FOR  A  NUMBER  OF  POLYURETHANE  FOAMS 


♦Taken  from  Reference  4 
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FIG.  6 


TYPICAL  STRESS- STRAIN  CURVES  FOR  3.5  LBS./FT : 
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FIG.  11  DENSITY  DEPENDENCE  OF  PLATEAU  STRESS  OF  RIGID  POLYURETHANE  FOAM 
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FIG.  12  STRAIN  RATE  DEPENDENCE  OF  PLATEAU  STRESS  OF  RIGID 
POLYURETHANE  FOAM 
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FIG.  13  VARYING  AREA  MITIGATOR 
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FIG.  15  IMPACT  CONFIGURATIONS 
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FIG.  16  TYPICAL  DATA  RECORDS  OF  STRAIN  IN  ROD 
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17  TYPICAL  THEORETICAL  FORCE-TIME  CURVES  FOR  CONSTANT  AREA  MITIGATORS 
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FIG.  18  TYPICAL  THEORETICAL  FORCE-TIME  CURVES  FOR  VARYING  AREA  MITIGATORS 
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APPENDIX  A 


NUMERICAL  INTEGRATION 


Equations  (A-l)  and  (A-2) 


ET  t 

J°T  “/  /  I 

*o  *'o 


dtdt 


_F _ 

(pgA) 


(A-l) 


1 

L 


T 

-I’ 


_F _ 

(pgA) 


(A-2) 


are  numerically  integrated  using  Equation  (A- 3) 


F  =  F  (e ,  t,  A) 


(A-3) 


as  a  constraint.  The  integration  is  performed  by  applying  a 
successive  approximation  method  to  a  simple  quadrature  formula, 
the  trapezoidal  rule,  Figure  A-l.  Equations  A-l  and  A-2  become 


£n  -  Z  Vn  -  5S  j1<Fi  +  pi-l>  (ti-ti-i) 


(A-4) 


[**  -  *1-1]  - 


(Pi  +  pi-i>  (ti  -  W 


A-l 
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for  any  time  (K  tn<T.  At  each  new  increment  in  time  tn  =  tn_^  +  A  t, 

an  estimate  is  made  for  the  mitigator  force  Fn.  Since  the  force  is 
not  changing  very  rapidly  with  time  for  much  of  the  calculation ,  this 
assumed  force  is  merely  the  Fn  at  the  previous  time  tn  , .  The  strain 
en  and  strain-rate  e  at  time11  t  are  determined  by  evaluating 
Equations  (A-4)  and  (A-5)  using  this  FR.  A  new  force  Fn  is  then 
calculated  from  Equation  (A-3)  using  the  new  e  and  e  .  nThe  iteration 
is  continued  until  the  change  in  Fn  becomes  liss  than  .1  percent  of 
Fjj.  The  results  at  time  tn  are  recorded,  the  time  is  advanced  one 
time  step  and  the  entire  procedure  is  repeated  until  F  drops  to 
zero. 
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APPENDIX  B 

COMPUTER  PROGRAM  LISTING 


The  computer  program  is  written  to  handle  the  case  of  a 
rigid  projectile  impacting  a  semi-infinite  elastic  aluminum  rod 
protected  by  a  rigid  polyurethane  foam  shock  mitigator.  It 
reduces  the  raw  experimental  data,  calculates  the  equivalent 
theoretical  results,  and  prints  and  plots  the  experimental  and 
theoretical  results. 
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PROGRAM  CBBOSCR ( I NPUT , OUTPUT , TAPE5= INPUT ,TAPE6=0UTPUT ,TAPE99 ) 

1  READ(5,2)  NTEST *KOUT tNSPEC 

2  FORMAT (315) 

NTEST  =  TEST  IDENTIFICATION  NUMBER 
NTEST  .LT.  0  —  STOPS  RUN 
KOUT  =  TYPE  OF  RUN 

KOUT  =1  —  EXPERIMENTAL 

KOUT  =  2  —  THEORETICAL 

KOUT  =  3  —  EXPERIMENTAL  AND  THEORETICAL 

NSPEC  =  SPECIMEN  IDENTIFICATION  NUMBER 
iFtNTEST.LT. 0)  GO  TO  9999 

READ! 5*22)  SPECL, SPECDl ,SPECD2*SPECR0, PROJWT* VELIMP, VELRB 
22  FORMAT ( 7F10»5 ) 

SPECL  =  SPECIMEN  LENGTH  (IN) 

SPECD1  =  SPECIMEN  DIAMETER  (IMPACT  FACE)  (IN) 

SPECD2  =  SPECIMEN  DIAMETER  (REAR  FACE)  (IN) 

SPECRO  =  SPECIMEN  DENSITY  (LBS/FT3). 

PROJWT  =  PROJECTILE  WEIGHT  (LBS) 

VELIMP  =  IMPACT  VELOCITY  { IN/SEC) 

VELRB  =  REBOUND  VELOCITY  (IN/SEC) 

PROJEO  =  ,5*PR0JWT*VELIMP**2/386» 

PROJEF  =  • 5*PROJWT*VELRB**2/386. 

CALL  CALCM1  ( 1 ,0. > 0. ,0 ,0. ,0.001 *0. *200. *5. ,5. , 

123HTRANSMITTED  FORCE  PULSE, -23» 

210HTIME  (SEC)  ,10,I1HFORCE  (LBS) ,11 ,0. *18 ) 

WR  I TE ( 6 ,3 )  NTEST 

3  FORMAT( 1H1,11HTEST  NUMBERiI7) 

WR ITE (6,4)  NSPEC, SPECL, SPECRO 

4  FORMAT ( 1H0 » 1 1HSPECIMEN  NO, 1 5/ 11HLENGTH  (IN),F10.3/ 

122HDENSITY  (LBS/CUBIC  FT),F10.1) 

WRITE ( 6 ,5 )  PROJWT , VELIMP ,PROJEO»PROJEF 

5  FORMAT ( 1H0  »10HPROJECTILE/12HWEIGHT  (LBS) ,F10.6/24HIMPACT  VELOCITY 
1( IN/SEC) ,F10*0/24H INCIDENT  ENERGY  ( IN-LBS) ,E10.4/23HREBOUND  ENERGY 
2  ( IN-LBS) ,E10. 4) 

PROJWT  =  PROJWT/386. 

GO  TO  ( 1000,2000,1000) , KOUT 
1000  CALL  EXPMT(PROJEO, PROJEF) 

IF(KOUT.EQ.I)  GO  TO  3000 

2000  CALL  THEO( VELIMP* PROJWT ♦ PROJEO, SPECL, SPECD1 *SPECD2  * SPECRO) 

3000  GO  TO  1 

9999  CALL  CALCM1(0,0.) 

STOP 

END 


no  n  r»  n 
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SUBROUTINE  EXPMT ( PROJEO ♦PROJEF ) 

DIMENSION  TIMEM ( 200 ) » FORCEM ( 200 ) 

READ ( 5  » 1 )  NPTS,CNVRT,CNVRF 

1  FORMAT (1I5*2F10«5) 

NPTS  =  NUMBER  OF  EXPERIMENTAL  POINTS 

CNVRT  =  CONVERSION  FACTOR  FOR  TIME  (INPUT  UNITS  TO  SECONDS) 

CNVRF  =  CONVERSION  FACTOR  FOR  FORCE  (INPUT  UNITS  TO  LBS) 

DO  100  I =1 »NPTS»4 

READ( 5*2)  TIMEM ( I ) , FORCEM ( I ) *TIMEM( 1+1 ) *FORCEM( 1+1) ,TIMEM( 1+2) , 
1F0RCEM( 1+2) *TIMEM( 1+3) »FORCEM( 1+3) 

2  FORMAT (12X»8F6*0) 

100  CONTINUE 

TIMEM( I )  =  DATA 
FORCEM ( I )  =  DATA 
FINT  «  0. 

DO  200  1=1 *NPTS 

TIMEM ( I )  =  CNVRT*TIMEM( I ) 

FORCEM ( I )  =  CNVRF*FORCEM{ I ) 

IF  (I.EQ.l)  GO  TO  200 

FINT  =  FINT+ ( FORCEM( I )**2+FORCEM( 1-1 )**2)*(TIMEM( I ) -TIMEM ( 1-1) )/2. 
200  CONTINUE 

ENGYTR  =  • 102*F INT 

ENGYAB  =  PROJEO-PROJEF-ENGYTR 

CALL  CALCM1(-NPTS»TIMEM, FORCEM, 0) 

WR I TE ( 6 ,3 ) 

3  FORMAT ( 1H-»20HEXPER I MENTAL  RESULTS) 

WR I TE ( 6 ,4) 

4  FORMAT ( 1H  ,10HTIME  (SEC)  , 10X , 11HFORCE  (LBS)) 

DO  300  1=1 »NPTS 

WR I TE ( 6 * 5 )  TIMEM(I) ,F0RCEM(I) 

5  FORMAT ( 1H  ,2(1PE13.3) ) 

300  CONTINUE 

WRITE (6, 6)  ENGYAB »ENGYTR 

6  FORMAT ( 1H0»24HENERGY  ABSORBED  ( IN-LBS ) »E10. 4/ 

127HENERGY  TRANSMITTED  ( IN-LBS ) ,E10.4 ) 

RETURN 

END 
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SUBROUTINE  THEO(UO*M1,PROJEO,L»D1*D2»RO) 

DIMENSION  TIMEC ( 200 ) *FORCEC (200) *STRANC(200) *RATEC( 200) 

2  *VPROJ ( 200 ) »VBaR(200) *ENGYAB(200) »ENGYTR(200) 

REAL  I *  1 0  *  1 1 » 1 2  »M1 *L 
READ (5*1)  NCASE, STPSIZ 
1  FORMAT ( I5,E15.5) 

NCASE  =  TYPE  OF  MATERIAL  MODEL 
NCASE  =  1  —  RIGID,  PERFECTLY  PLASTIC 

NCASE  =  2  --  RIGID,  PERFECTLY  PLASTIC,  STRAIN  RATE  DEPENDENT 

NCASE  =  3  —  ELASTIC,  PERFECTLY  PLASTIC 

NCASE  '  4  —  ELASTIC,  PERFECTLY  PLASTIC,  STRAIN  RATE  DEPENDENT 

STPSIZ  =  TIME  STEP  (SECONDS) 

J  =  0 
NPT5C  =  0 
Z  *  10.2 
D  =  STPSIZ 
T  =  D 
F  *  0. 

FI  =  0. 

F5  =  0. 

F9  =  0. 

E9  =  0. 

I  =  0. 

10  =  0. 

110  DO  160  N=1 *200 

II  =  1+ ( F+Fl ) *D/2 • 

12  =  I0+(I+Il)*D/2. 

R1  =  ( U0-F1/Z ) /L-I l/( L#M1 ) 

El  =  U0*T/L-Il/(L*Z)-I2/( L*Ml ) 

I F ( NCASE.GT .2 )  GO  TO  153 

CALL  R IGID ( NCASE, UO,L,D1,D2*RO, El *R1*F1) 

IF(Rl.LT.O.)  GO  TO  5000 
GO  TO  155 

153  CALL  ELASTC(NCASE,J,U0,L,D1,D2,R0,E1,E9,R1*S*F1»F9) 

IFtFl.LT.O.)  GO  TO  5000 
155  IF(ABS(1-F5/F1).LT..001)  GO  TO  170 
F5  =  FI 
160  CONTINUE 
170  NPTSC  =  NPTSC+1 
TIMEC(NPTSC)  =  T 
FORCEC (NPTSC )  =  FI 
STRANCt NPTSC)  =  El 
RATEC ( NPTSC )  =  R1 
VPROJ ( NPTSC )  =  U0-I1/M1 
VBAR ( NPTSC )  =  Fl/Z 
IF  (NPTSC. NE.l)  GO  TO  175 
ENGYTR(l)  =  • 102* FORCEC ( 1 ) *#2*D/2 • 

GO  TO  176 

175  ENGYTR(NPTSC)  =  ENGYTR ( NPTSC-1 )+. 102* ( FORCEC ( NPTSC ) **2 
1+ FORCEC (NPTSC-1 )**2)*D/2. 

176  ENGYAB(NPTSC)  =  PROJEO-. 5*Ml*VPR0J ( NPTSC ) **2/386. -ENGYTR (NPTSC ) 
IF  (NPTSC. EQ. 200)  GO  TO  5000 
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180  T  =  T+D 
F  =  FI 
I  =  II 
10  =  12 

IF(E1.LE.E9)  GO  TO  240 
E9  =  El 
F9  =  FI 
240  GO  TO  110 

5000  CALL  CALCM1 ( -NPTSC »T IMEC ♦ FORCEC  *0 ) 

WRITE  (6*2)  NCASE *S 

2  FORMAT  (1H-,19HTHE0RETICAL  RESULTS , 10X »4HCASE » 1 5 * 10X , 3HE  =»F10.0) 
WRITE  (6*3) 

3  FORMAT  (1H  ,10HTIME  (SEC)  .  10X ,  1 1HF0RCE  (LBS)*3X» 

214HSTRAIN  ( I N/ I N ) ,3X , 12HRATE  ( 1/SEC) »3X,14HVPR0J  (IN/SEC),3X* 
313HVBAR  ( IN/SEC) »3X»11HEA  ( IN-LBS) »3X ,11HET  (IN-LBS)) 

DO  100  K=l, NPTSC 

WRITE  (6*4)  TIMEC(K) ,FORCEC(K) ,STRANC(K) ,RATEC(K) ,VPROJ(K) , 
2VBAR(K) ,ENGYAB(K) »ENGYTR(  K) 

4  FORMAT  ( 1H  *8( 1PE15.3) ) 

100  CONTINUE 

RETURN 

END 
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SUBROUTINE  RIGID  (NCASE»U0  »L»D1 »D2»RO,El »R1 »F1 ) 
C  RIGID,  PERFECTLY  PLASTIC  MATERIAL  MODEL 

REAL  L 

CALL  AREA  (A,D1,D2»E1) 

IF  (NCASE.EQ.2)  GO  TO  550 
R5=U0/L 
GO  TO  560 
550  R5=R1 

560  CALL  STRATE(R0,R5, STRESS) 

FI  =  A*STRc33 

RETURN 

END 


SUBROUTINE  ELASTC(NCASE»J,U0,L,Dl»D2»RO,EltE9»Rl*S,Fl,F9) 
C  ELASTIC,  PERFECTLY  PLASTIC  MATERIAL  MODEL 
REAL  L 

CALL  AREA  (A,D1,D2,E1) 

CALL  MOD(RO,S) 

IF  (NCASE.EQ.4)  GO  TO  625 
R5=U0/L 
GO  TO  630 
625  R5=R1' 

630  IF  (J.GT.O)  GO  TO  650 

CALL  STRATE(RO,R5, STRESS) 

EYIELD  =  STRESS/S 
IF  (El. GT. EYIELD)  GO  TO  650 
FI  =  S*E1*A 
GO  TO  685 

650  IF  (E1.LT.E9)  GO  TO  675 
IF  (E1.GT..5)  GO  TO  675 
J  *  J+l 

IF  (Rl.LT.O.)  GO  TO  675 
CALL  STRATE(RO,R5, STRESS) 

FI  *  A*STRESS 
GO  TO  685 

675  FI  =  F9-S* ( E9-E1 ) *A 
685  RETURN 
END 
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SUBROUTINE  MOD(RO,S) 

DENSITY  DEPENDENCE  OF  MODULUS 

S  =  230.*RO**1.7 

RETURN 

END 


SUBROUTINE  STRATE ( RO,R5 , STRESS ) 

STRAIN  RATE  AND  DENSITY  DEPENDENCE  OF  STRESS 
STRESS  =  (8.69+.306*ALOG(R5) )*R0**1.69 
RETURN 
END 


SUBROUTINE  AREA  (A,D1,D2»E1) 
STRESS  AREA 

A=»7854#(Dl+(D2-Dl)*El/»5) *#2 • 

RETURN 

END 


